Babinet's principle in the Fresnel regime studied using ultrasound 
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The diffraction of ultrasound by a circular disk and an aperture of the same size have 
been investigated as a demonstration of Babinet's principle in the Fresnel regime. The 
amplitude and the phase of diffracted ultrasonic waves have been measured and a 
graphical treatment of the results is performed. The results verify Babinet's principle. It 
is also found that the incident wave is indeed nil behind the phase of the wave passing 
through on the central axis of a circular aperture. This paradox has previously been 
regarded as a defect of Fresnel 's theory. 
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I. INTRODUCTION 



The high intensity spot just behind a circular opaque obstacle, known as the 
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Poisson's bright spot (the spot of Arago), ' arouses curiosity and can be used to gain the 
attention and interest not only of students of physics but also of those of general science 
classes. The additional observation of the complementary circular aperture provides 
materials for quantitative discussions on the mathematical theory of waves. 

A general description of Babinet's principle states that the sum of the diffraction 
fields U\ and U2 behind two complementary objects is the field Uq observed without the 
diffracting objects. " 

U Q = U l + U 2 . (1) 
This principle is generally fulfilled within the framework of the scalar theory of 
diffraction. Verification of the Babinet's principle has been almost exclusively 
performed with light. The diffracted intensity patterns (7 = UU*) are compared in the 
Fraunhofer regime 4 ' 7 (source and observation plane situated very far from the diffracting 
object). In the Fresnel regime (source and observation plane at a finite distance from the 
diffracting object), the intensity patterns are different from each other and are compared 
through calculation. This is mainly because in optics the wavelength is too short and 
the frequency too high to observe the phase. In the microwave region, the wavelength 
becomes a few cm. It is possible to make use of this advantage, e.g. the phase 
dependence of the dielectric transient absorption allows observation of a phenomenon 
several orders of magnitude faster than the time resolution which is determined by the Q 
value of the cavity. 9 The wavelength of ultrasound is also several mm and the period is a 
few tens of microseconds. Then, it is possible to manipulate and observe the phase as 
well as the amplitude. This enables Babinet's principle to be studied directly by 
comparing the amplitude and the phase. With ultrasound measurements the essence of 
the theory can be seen by simply drawing vectors (phasors) in the complex plane, 
without complicated calculation. 

We will observe diffraction by a circular aperture having a diffracted field U a and a 
complementary circular disk and study Babinet's principle, Uq = U a + by using 
ultrasound. The apparatus is inexpensive and can be assembled by hand. The 
measurement is as simple as learning how to use the oscilloscope. A graphical treatment, 
the vibration spiral, based on Huygens-Fresnel diffraction theory, is used to analyze the 
observed amplitude and phase. 

A curious feature of both the graphical approaches and mathematical analysis is a 
n/2 advance of the phase for the wave passing through on the central axis of a circular 
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aperture with respect to the incident wave Uq. This phase difference has been regarded 
as a defect of Huygens-Fresnel diffraction theory. However, the phase difference has 
never been actually measured as far as we know. We will measure the phase due to a 
small aperture with a diameter close to the wavelength being used. The vibration spiral 
for the circular aperture will be investigated in detail. 



II. THEORETICAL BACKGROUND 

Huygens' principle explains the propagation of a wave as follows: every point of a 
wave-front may be regarded as a source of secondary spherical wavelets, and the 
wave-front at any later instant is constructed as the envelope of these wavelets as shown 
in Fig. I. 5 ' 10 Fresnel constructed a diffraction theory based on Huygens' principle 
together with the principle of interference. We will discuss the theoretical basis within 
two alternative approaches of a graphical and mathematical treatments. 

A. Graphical treatment of diffraction by a circular object 

A graphical treatment, based on Huygens-Fresnel diffraction theory, gives a clear 
physical picture of the origin of the diffraction pattern. 4 ' 6 ' 7 Consider a point source S 
located at a distance zq from a circular aperture and the diffracted field U a at point P on 
the axis at a distance z from the aperture. Consider now a set of imagined spheres 
centered at P of radii z + A/2, z + 2 A/2, z + 3A/2 and so on, differing by A/2 as shown in 
Fig. 1 . Their intersections at the aperture with a spherical wave front emanating from 
the source S form circular zones, F\, F2, F3 and so on. The zones are called Fresnel or 
half-period zones. There is a n phase difference between consecutive zones, and 
therefore, the contributions of the successive zones to the diffracted field are alternately 
positive and negative. Figure 2(a) shows the amplitude diagram in the complex plane 
when the first half-period zone is divided into eight subzones. 4 The phasor a\ represents 
the contribution from the first subzone. The phasors a-i and due to the second and 
third subzones are added. The phase difference 8 between each successive zone will be 
7t/8. Adding all eight subzones gives the phasor OG as the resultant amplitude from the 
first half-period zone. Similarly, the contribution from the second half-period zone is 
GH, and OH is the sum of the first two zones. Successive half-period zones give the rest 
of Fig. 2(a). By increasing the number of subzones indefinitely, the vibration spiral is 
obtained as shown in Fig. 2(b). U a is represented by phasor OD. The spiral eventually 
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approaches C; the sum of all the half-period zones covers the whole spherical wave, 
which should be the phasor OC representing the unobstructed field Uq. 

In Babinet' principle, Uo=U a + Ud, where Ud is DC. The spiral is practically 
circular for many turns. Huygens-Fresnel diffraction theory predicts that the Poisson's 
bright spot appears behind a circular disc, its amplitude being almost equal to Uq. The 
argument here also applies for the plane wave geometry. 

The spiral also allows one to determine directly the resultant amplitude and the 
phase due to any fractional number of zones. The phase can be estimated by simply 
calculating As, the difference between the path length along the axis and that via the 
edge of the circular object. The angle cp (ZOCD) on the spiral (Fig. 2) is given by 2k X 
Ab IX. Since Ab can be larger than X, cp can exceed 2k. 

Figure 2 applies for the time t = 0. The phasors revolve clockwise with the same 
angular velocity co. The triangle representing Uq, U a and Ud revolves as a rigid frame. 
The simple harmonic wave is presented as the projection on the coordinate of the real 
part. We adopt here the custom in optics of representing an advance of phase by a 
clockwise rotation of the amplitude phasor. 

A problem arises. The resultant phasor Uq turns out to be k/2 behind the wave from 
the center of the zone system, whereas, according with Huygen's principle, a wavelet 
emitted at O (Fig. 1) in phase with the incident wave Uq should arrive at P still in phase 
with Uq. This paradox has been regarded as a defect of Fresnel's theory resulting from 
the approximations made therein. 4 ' 7 

B. Mathematical treatment in relation to the vibration spiral 

We will discuss now the scalar theory of diffraction in relation to the vibration 
spiral. Kirchhoff gave a mathematical expression for Fresnel's wave theory on the basis 
of Green's theorem. Fresnel-Kirchhoff diffraction theory provides an accurate solution 
in many instances and used most widely. The diffraction patterns for the circular 
aperture and disk are given by Lommel functions, 1 which are described in terms of the 
Bessel functions. We consider for simplicity the diffraction on the axis. 

Fresnel-Kirchhoff diffraction for a point source and a circular aperture in polar 
coordinates can be expressed in terms of the incident wave Uq(P). One obtains for point 
P on the axis, 6 




(2) 



U (P) = A exp {- i(o)t -kr- 0)}l r , 



(3) 
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where k = 2n/X, Ao is the amplitude of the incident wave, 90 is the value of the phase 
difference at P between a wave that comes from the edge of the aperture and that which 
follows the z-axis. The inclination factor Q(ty) = '/i(cosa + cosP) is a smooth 
monotonically decreasing function, a and P are angles that the diffracted and incident 
waves, respectively, would make with a normal unit vector n perpendicular to the 
diffraction plane (Fig. 1). The -i appearing in Eq. (2) accounts for a n/2 phase shift of 
the diffracted wave relative to the incident wave. 5 " 7 ' 10 The phasor integral 

l{cp,) = ^Q{<P)e iv dcp (4) 

gives the vibration spiral. For a few turns Q(ty) is close to 1 . Therefore, the spiral may 
be regarded as a circle, /(cpo) = l+i, 2i, -l+i and for the phase difference cpo = n/2, n, 
3n/2 and 2n, respectively, and the points for /(cpo), multiplied by^o, will lie on the circle 
in Fig. 2(b). In the limit of infinite aperture radius, /(cpo) becomes i, and U a (P) 
approaches Uo(P), the incident wave. The resultant phase from Eq. (2) is -n/2 for a 
wave that follows the z axis and for no aperture. The vibration spiral in Fig. 2 should 
be turned n/2 clockwise, then Uq is on the real axis. 

However, Fresnel-Kirchhoff diffraction presents inconsistencies in the boundary 
conditions and is valid when the diffraction object is large compared with the 
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wavelength (D»X) but small compared with zo and z (zo, z » D). ' 

Rayleigh-Sommerfeld diffraction theory gives a rigorous expression for the scalar 
field. 11 This theory has not been discussed in relation to the vibration spiral. Closed 
solutions for the diffracted field for the spherical wave by a circular aperture and a disk 
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along the axis is given as, 



exp{//V(z n + z)} z exp{ik(Z n + Z)\ 

—, \ + exp(*;r) ; " 

(z +z) FV 'z (Z 0+ Z) 



(5) 



zexp^ + Z)} 

z (z + z) 

where Zo is the distance between the edge and the point source on the axis and Z is the 
distance between the edge and the observation point on the axis (Fig. 3). These simple 
formulas are valid anywhere along the axis behind the object. The first term in Eq. (5) is 
the incident spherical wave Uo(z) and the second term is -Ud(z), ( exp(/;r) = -1 ). We 
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have the Babinet's principle, Uo= U a + on the axis. 

Now we move into the Fresnel regime (zo, z > D) and look at Ud{z) as a function of 
D at a fixed z on the axis. Eq. (6) is rewritten as, 
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tj (*\ a z z ° +z Q *p{ ik ( z o +z )l n\ 

U d ( z ) = A o- - ^ 7 7 - x • (7) 

Z Z + Z Jz + z J 

The amplitude for a disk is given by (z/ Z)-(z + z)/(Z +Z)- A and is almost equal 
to v4o when £) is not too large. The phase difference from the incident wave is given by cp 

= 2nA B /X = 27r{(Z + Z) - (z + z)}/ X . Let us put the endpoint of phasor U& onto that of 

Uo and see the trace of the starting point of which is also the end point of U a , on the 
complex plane by increasing D. Then we have the vibration spiral as shown in Fig. 2(b), 
except Uo is on the real axis. For a very small aperture, Ab is small and Ud becomes 
close to Uo. Therefore, the resultant phasor U a = Uo - Ud has a phase n/2 ahead from Uo. 
Fresnel zones can be obtained as (z + z) + A B = (z + z) + n ■ A 1 2 , where n=l ,2, 3 

and so on. 

The diffracted field for the plane wave by a circular aperture along the axis is 
expressed simply, 13 ' 14 

U a (z) = A z\ eMikz) - eMikZ) ]. (8) 
z Z 
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An one-dimensional integral over the rim of the aperture and a conventional 
two-dimensional integration over the area of the aperture 14 yield the same result. 

As we have seen, the theories all show a n/2 phase difference between the incident 
wave and the wave passing through on the central axis of a circular aperture. 



III. EXPERIMENTAL METHOD 

The experimental setup using ultrasound is shown schematically in Fig. 3. The 
transmitter T (Nicera T4016) was driven at 40 kHz by an oscillator. This gives a 
wavelength X of 8.58 mm at room temperature. The frequency was monitored by a 
frequency counter. The transmitter and the center of the diffracting objects were placed 
on the axis so that the ultrasonic path has cylindrical geometry about the normal to the 
object through its center. The receiver R (Nicera R4010A1) was chosen because it has a 
wide signal accepting angle (95 ° FWHM) and was set on a rail placed perpendicular to 
the axis and scans the radial distribution of the signal beyond the circular aperture or a 
disk. The distances zo and z are from the diffracting object (the aperture or the disk) to 
the source and the scanning line respectively. A transmitter with narrow emitting angle 
(50 °FWHM) was chosen in order to avoid undesirable reflection mainly from the table. 
The distance zo was set to be larger than z in order to have a uniform sound pressure at 
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the aperture plane. The signal from the receiver was put into chl of an oscilloscope and 
the voltage to the transmitter was monitored in ch2. The voltage used was typically 2-4 
Vp.p. Then the amplitude of the received signal and the phase difference with respect to 
the signal in ch2 were observed. The amplitude and the phase were registered as a 
function of the radial distance r. The diameters of the transducers were 7 mm and 6 mm 
for T4016 and R4010A1, respectively. 

A. Circular aperture and disk for study of Babinet's principle 

A circular aperture and a complementary disk of 10 A, in diameter were prepared 
for the study of Babinet's principle. A circular aperture of 86 mm in diameter cut into 
0.6 mm thick A4 size hard paper was placed on a cork board of 900 mm wide, 600 mm 
high and 6 mm thick with a circular hole of 1 1 8 mm in diameter. The cork board acts as 
a frame supporting the aperture and also absorbs the sound. A circular disk of the same 
diameter was also prepared. A circular disk taken from the same paper was backed by a 
circular cork disk of 69 mm in diameter and 4 mm thick. The disk was supported by a 
carbon rod of 4 mm in diameter. The apparatus was set up on a 5 cm-square-ruled mat 
which covered the desk to guide alignment. The heights of the centers of the transmitter, 
the receiver, the aperture and the disk were set to be 340 mm. 

The receiver was put on a rail placed 858 mm (100 X) away from the transmitter so 
that the radial distribution of the diffracted wave can be obtained. The position on the 
vibration spiral was selected by altering the position z of the aperture as shown in Table 
I. 

B. Small aperture for phase study of vibration spiral 

According to the vibration spiral, the phase difference A9 observed on the axis for 
a circular aperture field U a relative to the unobstructed field Uq will be -n/2 in the 
zero diameter limit, and -n/4, 0, and n/2, respectively, for an aperture containing a 
quarter-period zone, a half-period zone, and one period zone (Fig. 2). Given the 
relatively long wavelength of the ultrasound, the diameter of the first Fresnel 
half-period zone can be several cm. Therefore it is possible to observe effects within a 
small fraction of the half-period zone. A small aperture, which will contain only a tiny 
fraction of half-period zones, was prepared. The phase differences comparing with and 
without the aperture were measured. 

We used a small circular aperture of 1 AX in diameter and used the graphical 
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treatment. A4 size hard paper of 0.6 mm thick with a 12mm diameter hole at the center 
was placed on the aperture system used in experiment A. The zo + z distance was also 
set to be 100A.. The fraction of half-period zone covered was changed by altering the 
position of the small aperture. The aperture was set between z = 17X- 38A,.The fraction 
can be as small as 1/30 - 1/50 half-zone. The phase and the amplitude of the 
unobstructed wave U$ on the axis were compared with the corresponding phase and 
amplitude of the wave U sa with the small aperture. The phase was also observed off the 
axis as a function of radial distance r at z = 17A, and 38A. The position cp on the vibration 
spiral and the expected phase difference AG (= -90°+ cp/2) on the axis are shown in Table 
II. Sublabels 1, 2 and 3 stand for the distances z of 17X, 26. 3A, and 38A, respectively. 

IV. RESULTS 

A. Babinet's principle 

The relationship between Uo, U a and Ud for the 10A, circular aperture and the disk 
placed 17X, 26.3X and 38A, from the receiver set on the axis are shown in Fig. 4. The 
distance between the transmitter and the receiver was 100A,. The experimental 
parameters are shown in Table I. The end points for phasors U a and Uo - Ud are plotted 
as open and closed symbols. Measurements of U a and U& can be compared with theory 
independent of each other in this plot. The phase 6o for no diffracting object was taken 
to be n/2 following the custom of representing the Fresnel vibration spiral. If the 
assumption in graphical treatment is true, the measured points lie on the spiral. The 
amplitude values for the aperture A a and the disk Ad are normalized to Aq. The cp values 
calculated by the path difference Ab for each setup are shown with bars on the circle. 
The aperture contains less than two half-period zones with these setups. The scattered 
points are results for different measurements. The arrows show the average values forz 
= 17A,. The main contribution to the error comes from errors in the phase. The amplitude 
for Ud is observed to be larger than that for Uo and points are mostly outside the circle 
except for z = 17X. On the other hand, most points for U a lay inside of the circle. The 
reason may be due to the finite thickness of the edges, and the limited radial width of 
incident ultrasound. Although, there still are small discrepancies in points for U a and Ud, 
the agreement between the graphical treatment and experimental results is good. It has 
also been shown that the point moves counter clockwise on the spiral as A B increases 
when the object distance is altered from the center to the receiver. 
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The A a values on the axis can be between and twice Aq depending on the position 
z of the aperture. In the present setups, two A a maxima were observed at around z = 44X 
and 56A, where A B is the same. A a was observed to be about 1.9 Aq and Q a was almost 
equal to 60 at the maximum. An A a minimum was observed at z = \4X. Q a changed 
abruptly from being behind to ahead of 9o as the aperture approaches the receiver in the 
vicinity of the minimum. 

The amplitude and the phase measurements were also performed for off-axis 
positions. The typical radial diffraction patterns, the amplitude A and the phase 9, 
obtained for Uq, U a and U& are shown in Fig. 5. There are five characteristic positions 
shown in Fig. 5. One is the axis, and two are geometrical edges (vertical broken lines) 
which lay symmetric about the axis. Other positions M\, where all phases come close, 
correspond to positions for the eclipse of the first half-period zone, (dot-dash lines). 
When the observation point P moves out from the center and goes into the geometrical 
shadow of the aperture, the Fresnel zone system follows. As P moves further, the total 
eclipse of a half-period zone occurs. 

Phases also come close to each other at the geometrical edges. The position of the 
geometrical edge and the position where three phases coincide differ a little from each 
other. The difference may be due to the finite sizes of the transducers. 

Aq and Go change smoothly as expected over the whole range observed. Observed 
Aq values are considerably flat and almost constant. 60 lags as the receiver moves out 
from the axis and the distance from the source increases. The parabolic curve shows the 
phase shift of Uq, 27t{TR'-(z q + z)}/ A , calculated by the difference in 

transmitter-receiver distance between the receiver on the axis (zq+z) and on a radial 
position (77?'). The curve represents the phase of a spherical wave originated at the 
transmitter observed on a line perpendicular to the axis as a function of radius r. 
Observed 9o agree well with the calculated curve as expected. 

The amplitude A a for the aperture increases slowly when the receiver approaches 
the geometrical edge from outside. Two A a maxima lay symmetric about the axis inside 
the geometrical edges. A a on the axis depends on z as mentioned above. The amplitude 
A^ for the disk is almost constant and the same as Aq, outside the shadow except a small 
minimum near Mi, until the receiver approaches the geometrical edge. Then it starts 
decreasing as the receiver comes close to the geometrical edge. Structures are seen in 
the shadow. Poisson's bright spot appears on the axis surrounded by small maxima. The 
amplitude of the spot is almost equal to Aq as predicted by the graphical treatment. 

Change in 6 a is slow inside the edge, except in the vicinity of the center. However, 
the phase Q a changes rapidly outside the edge (in the geometrical shadow). Outside the 
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edge, Q a goes behind as r increases until it approaches position M\, then Q a starts the 
reversal movement. The phase Q a goes ahead instead of behind. After passing through 
Mi, Q a again goes behind as r increases. However, A a is small and the changes inA a are 
not evident in this region. The same behavior is observed at the boundaries of the 
successive outer zones. 

The phase 6d is almost the same as 60 outside the geometrical shadow. Inside the 
shadow, 9d goes behind as the receiver goes inside. The 6d changes rapidly in the 
vicinity of Ad minima. 6d on the axis differs by about 2n from 60 with this setup. One 
cannot tell the difference by only observing phases for 9d and 9o near the axis. However, 
the difference is clear by observing over a wide range of r. This is also seen in the 
vibration spiral as shown in Fig. 4 where U& has rotated about 2n. 

The amplitude patterns for aperture A a and disk A^ are very different from each 
other as shown in Fig. 5 and one cannot tell the relation of U a and by just looking at 
those patterns. However, Babinet's principle is satisfied when one compares the 
amplitudes and phase, off the axis. The phasors representing Uq, U a and U& off the axis 
do not lie on the circle. Still they make a triangle as expected by Babinet's principle. 
However, it was found that the phases tend to coincide or to have a n difference at 
'characteristic' points. Three phases coincide just outside of the geometrical edge and^o 
is given by A a + A&. The same relation is also observed at Mi, although the error is large 
because A a is very small. At the maxima of A a , 6d differs by n from 9o and Q a and the 
patterns show A =A a - A d . 

B. Phase difference 

The amplitudes and phases obtained for the 12 mm circular aperture are plotted on 
a vibration spiral as shown in Fig. 6. The end points for phasor U sa are plotted with 
different symbols according to z. Arrow U sa shows the average value for all z positions 
and Usa3 shows that for z = 3SX. The phase difference 9 between the path via the center 
of the small aperture and via the edge was very small as shown in Table II. Only a tiny 
fraction of a half-period zone appears in the small aperture. The 9 values calculated by 
the path difference Ab for z = 17A, and 3 8 A, are shown with bars on the arc. The observed 
9 agreed well with the calculated values. However, the observed end points for U sa .i and 
U sa 2 lay inside the circle as observed for U a in Experiment A; consequently, the values 
of |A6 sa | are observed a little smaller than the theoretical values. The reason may be due 
to the finite size of the receiver and errors in alignment. Influences of these are expected 
to be larger for smaller z. On the other hand the points for C/ sa 3 lay mostly on the arc and 
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the A9 value obtained was -85° and was close to the theoretical value of -88°. It is 
clearly seen that the wave from the center of the zone system is indeed nil ahead of the 
incident wave as expected from Huygens-Fresnel diffraction theory. 

The radial diffraction patterns for the amplitudes A sal and^ sa 3 obtained for a small 
aperture at z = 17A, and 38X, respectively, are shown in Fig. 7 together with^o- Vertical 
lines represent the geometrical edges. Abroad peak with a small structure at the center 
was observed for A sa \. The amplitude on the axis is almost 1/10 that observed for the 
unobstructed wave. The pattern for A sa i, is broader and the structure is not clear due to 
the smallness of the signal. The phases 9 sa i and 9 sa 3 observed at z = \1X and 38X, 
respectively, are also shown in Fig. 7. The dot-dash and broken parabolic curves are 
phase shift of U sa , IniOR'-z)/ A + Ad , calculated for the geometrical path difference 
between the center of the aperture to the receiver on the axis (z) and the center of the 
aperture to the receiver at r (OR'), with an added phase difference AG of -74° or -90°. A 
value of -74° was chosen as the average value for 6 sa observed on the axis and -90° is 
the value for the small aperture limit. The measured phase agrees well with either of the 
curves. The result also shows that the unobstructed wave Uo is indeed almost 90° behind 
the wave due to a very small aperture. The shape (width) of the parabolic curve is 
determined by the position of the origin, the center of the small aperture here, of the 
spherical wave. 



V. DISCUSSION 

The measurement of phase in addition to the amplitude makes it possible to verify 
Babinet's principle directly in the Fresnel regime. The observation of phase also shows 
some characteristics of diffraction otherwise difficult to notice. The zone system has 
been studied and various applications have been made such as zone lenses in the 
X-ray region. 7 ' 10 ' 15 The present observation of a rapid phase change at an eclipse of a 
half-period zone shows another aspect of diffraction. 

The phase and amplitude observed on the axis are simply predicted by the path 
difference between the wave along the axis and the wave via the edge of the aperture 
independent of the position of the aperture or the disk. The wave observed for a small 
aperture, which contains only a fraction of a half-period zone, has a phase always k/2 
ahead of the wave without the aperture, provided that the aperture is small enough and 
is not too close to the source or the observation point. 

The physical meaning of the origin of the n/2 phase difference is not fully 
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understood. The presence of the factor -i in Eq. (2) gives the necessary n/2 phase 
difference. However, it indicates that the secondary sources oscillate a quarter of period 
ahead of the primary disturbance. 5 ' 7 On the other hand, Huygens' principle states that 
the wavelets oscillate in phase with the primary wave. The parabolic curve 60 in Fig. 7 
represents a spherical wave centered at the transmitter, while 9 sa i or 9 S£I 3, after passing 
through the small aperture, represent spherical waves centered at the small aperture. The 
latter can be regarded as an envelope of Huygens' spherical wavelet originated at the 
small aperture but with the phase shifted ahead by n/2. 

A secondary wavelet approach by Huygens-Fresnel and an edge diffraction 
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approach by Young are two fundamental physical models of diffraction ' that could 
provide insight into the n/2 phase shift. The Huygens-Fresnel approach expresses the 
disturbance at P as the sum of all contributions from every point in the aperture as 
shown by the vibration spiral. The edge gives no contributions. An edge diffraction 
approach, the Maggi-Rubinowicz theory, describes the field behind the object as the 
superposition of the incident wave and the boundary diffraction wave (U B , wave 
originated at the edge shifted by n) in the illuminated region, and, in the geometrical 
shadow, just as the boundary diffraction wave, this is. ' 

U a (r) = A ex P^ r ) + jj b ( r ) when P is in the direct beam, (9) 
r 

U a (r) = +U B (r) when P is in the geometrical shadow. (10) 

Rubinowicz showed that the Kirchhoff diffraction integral (the mathematical expression 
of the Huygens-Fresnel approach) can be divided into the incident wave and the 
boundary wave, which is a line integral over the rim of an aperture, according to 
Young's approach. 3 ' 5 ' 16 The two approaches are considered as mathematically equivalent 
treatments. The edge diffraction theory, however has not been proved beyond the limits 
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of validity of the Kirchhoff theory, (zo, z» D, D»X). ' The physical reality of the 
boundary diffraction wave has been demonstrated by Ganci at the straight edge using 
the Young's double slit. 16 

We have tried to detect the boundary diffraction wave by observing the phase of 
the wave diffracted by a small aperture. The first terms in Eqs. (5) and (8) are the 
incident wave and the second terms (-Ud) can be regarded as the boundary diffraction 
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wave Ub. ' Then on the axis, we have U B = -Ud in Eqs. (5) and (8). As discussed in 
II.B., Uo and Ub have almost the same amplitude and a phase difference of n. Then, U a 
becomes n/2 ahead of Uo. The phase relation is also shown in the present measurements. 
When P moves deep into the shadow, Uo vanishes and only Ub remains. Then, U a 
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should have a phase difference of n in the shadow. However, the phase difference for U a 
remains n/2 both on the axis and in the geometrical shadow as shown by the parabolic 
curves 9 sa i and 6 Sfl 3 in Fig. 7. The reason may be that the "shadow" is not defined 
clearly when D becomes small. 7 The present experimental setup is beyond the validity 
of the Kirchhoff theory (zq, z» D »X). Sophisticated measurements like Ganci's 
experiment 16 may be needed to observe the physical reality of the boundary diffraction 
wave. 

The implemented apparatus was intended as a table-top instrument for the student 
laboratory and is far from ideal. However, it is still capable of showing fundamental 
aspects of the theory of diffraction. Ideally one needs a uniform wave and a huge room 
with no reflection in order to conduct a deeper investigation of the wave diffraction 
theory including the phase. In the classroom, however, satisfactory results can be 
obtained by setting the transmitter and the receiver close to each other with a proper 
alignment. Since Babinet's principle can be studied directly, the diffracting object need 
not be circular. An annular shape, 14 ' 15 an elliptical form or squares, in fact any form can 
be used. 



VI. CONCLUSIONS 

Babinet's principle in the extreme Fresnel regime has been directly verified by 
observing ultrasound waves diffracted by a circular aperture and a complementary disk. 
The phase as well as the amplitude were measured and analyzed by a graphical 
treatment. It has been also found that the wave without the diffracting objects is indeed 
n/2 behind the wave from the center of the zone system, which has been regarded as a 
defect of Fresnel's theory. The n/2 phase difference appears also in Fresnel-Kirchhoff 
diffraction, Rayleigh-Sommerfeld diffraction and even in the edge-diffracted approach 
by Young. 
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Table I. Experimental parameters for the study of Babinet's principle. Values for AJAq 
are obtained by closed solutions, Eqs. (5) and (6), for Rayleigh-Sommerfeld diffraction 
theory. 12 



frequency 

wavelength X at room temperature 
D 

Zq + Z 



40.0 kHz 
8.6 mm 
10 X 
100 X 



z 

A B 

A d /A 



(X) 
(X) 
(°) 



17 
0.87 
313 
0.95 



26.3 
0.64 
230 
0.98 



38 
0.53 
190 
0.99 



Table II. Experimental parameters for the study of the vibration spiral. AG is the phase 
difference from the incident wave. 



frequency 

wavelength X at room temperature 
D 

ZQ + Z 



40.0 kHz 
8.6 mm 
1.4 X 
100 X 



sublabels 



sal 



sa2 



sa3 



A B 

A9 (=-90+q>/2) 

A d /A 



(X) 
(X) 

(°) 
(°) 



17 

0.017 
6.3 
-87 
1.00 



26.3 

0.013 

4.6 

-88 

1.00 



38 

0.010 
3.7 
-88 
1.00 
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w 
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1 . Huygens' principle and construction of Fresnel's half-period zones (F\, F 2 , F3) 
on a spherical wave front at a circular aperture. A spherical wave front Wis 
originated at S. O is the center of the Fresnel zone and P is the observation point. 
Huygens' secondary wavelets originating on W combine to form a new wave 
front W. n is a normal unit vector perpendicular to the diffraction plane. 




2. Vibration spiral for Fresnel half-period zones of a circular aperture. 

(a) half-period zones are divided into 8 subzones, (b) Continuous subdivision of 
the zones. 
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Fig. 3. Schematic view of the experimental setup with a circular aperture. The data 
for the vibration spiral (Fig. 4) were taken on the axis (r = 0) by changing z. 



Im 




Fig. 4. The amplitude and phase relation obtained for unobstructed wave Uo, a 

circular aperture U a and a circular disk The receiver was set on the axis. The 
phase 6o for Uo was taken to be n/2. The end points for phasors U a and Uo - 
are plotted with open and closed symbols, respectively. The (p values calculated 
for each setup are shown with bars on the circle. 
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5. The diffraction patterns for the amplitude A (bottom) and the phase 6 (top) as 
a function of the radial distance r from the axis observed for complementary 
circular aperture and disk. The vertical broken lines (Eq) show the geometrical 
edges and dot-dash lines (Mi) represent positions for the eclipse of the first 
half-period zone (see the text). The phase range goes from below 0° up to 720° 
because of the continuity of the wave. 
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6. The amplitude and phase relation obtained for a small aperture of 1 AX in 
diameter. The receiver was set on the axis. The area enclosed by a square is 
displayed on an expanded scale on the right. The end points for phasors U sa \, 
U S a2, and U sa 3, at distance z of \7X, 26.3X and 38X, are plotted with triangles, 
squares and circles, respectively. Arrow U sa shows the average value for all z 
positions and U sa 3 shows that for z = 3SX. The (p values calculated for z = 17A, 
and 38A, are indicated with bars on the arc. 
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. 7. The diffraction patterns for the amplitude A (bottom) and the phase 9 (top) as a 
function of the radial distance observed for a small aperture of 1 AX in diameter. 
Sub labels 1 and 3 stand for the distances z of 17A, and respectively. 
Parabolic curves show phase differences calculated for the geometrical path 
difference between the center of the small aperture to the point on the axis (z) 
and the center of the aperture to the receiver (OR') with an added phase 
difference of -74° or -90° (see the text). The vertical lines (Eq S o) show the 
geometrical edges. 
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